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SUMMARY 

The  performance  of  arrays  of  receivers  can  be  improved  by 
optimisation  of  the  weighting  vectors  used  to  form  the 
output  beams . 

Three  optimisation  criteria  are  considered:  minimisation 
of  total  output  power,  least  mean  square  fit  to  a  signal 
vector,  and  maximisation  of  signal-to-noise  ratio.  For 
practical  purposes  it  is  often  necessary  to  apply  some  form 
of  constraint  to  the  weighting  vectors.  The  optimum 

narrowband  array  weight  vector  is  derived  for  each  of  the 
three  criteria  when  an  arbitrary  set  of  linear  constraints 
are  imposed  on  the  weight  vector. 

The  equality  of  the  derived  weight  vectors  in  most  cases  of 
practical  interest  is  proved. 

Finally,  some  examples  are  considered  which  show  how  a  number 
of  recently  proposed  weighting  vectors  can  be  derived  as 
special  cases  of  the  constrained  solutions. 
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1.  INTRODUCTION 

The  choice  of  array  shading  weights  modifies  the  performance  of  a  delay  and  sum 
beamformer.  In  conventional  beamforming  the  choice  of  a  set  of  weighting  vectors 
involves  a  tradeoff  between  the  width  of  the  main  beam  and  the  height  of  unwanted 
sidelobes. 

If  the  noise  field  is  known  *a  priori*  the  performance  can  be  optimised  in  some 
sense  by  a  correct  choice  of  the  weighting  vectors. 

This  choice  of  an  optirnmn  set  of  array  shading  weights  has  received  considerable 
attentionCref. 1,2,3)  and  the  weight  vectors  are  often  chosen  to  satisfy  optimisa¬ 
tion  criteria  such  as: 

(1)  minimising  the  total  output  power,  or 

(2)  a  least  mean  square  fit  to  the  signal  vector,  or 

(3)  maximisation  of  the  signal-to-noise  ratio 

In  general  a  single  point  constraint  of  a  fixed  (non  zero)  response  to  signals 
from  the  steered  direction  is  incorporated  and  under  this  constraint  the  three 
optimisation  criteria  result  in  the  same  weight  vector(ref . 1) .  This  constraint 

ensures  that  signal  output  power  increases  as  the  input  signal-to-noise  ratio 
increases .  A  further  extension  has  been  to  derive  optimum  weight  vectors  which 
satisfy  an  arbitrary  number  of  linear  constraints (ref . 2,3) .  An  example  of  this 

is  the  use  of  derivative  constraints  to  maintain  the  main  lobe  response(ref , 4)  and 
hence  prevent  cancellation  of  desired  signals  close  to  the  steered  direction. 

The  optimisation  criterion  used  in  these  cases  has  been  that  of  minimising  the 
total  output  power(ref ,2,3)  and  maximising  the  signal-to-noise  ratio(ref .4) . 

In  this  paper  the  optimum  weight  vectors  for  multiple  linear  constraints  are 
derived  for  the  three  different  optimisation  criteria  discussed  above.  Under 
some  very  mild  restrictions  on  the  choices  of  constraints,  it  is  shown  that  all 
three  criteria  result  in  the  same  weight  vector  for  a  given  set  of  constraints. 

The  frequency  domain  formulation  is  used  to  show  this  but  the  results  can  easily 
be  extended  into  the  time  domain. 

Some  special  examples  are  studied  to  show  how  the  general  solutions  reduce  to  a 
number  of  processors  that  have  appeared  in  the  literature(ref. 5,6,7) .  In 
particular  when  the  problem  is  overcorstrained  (ie  the  number  of  constraints  is 
greater  than  the  number  of  receivers)  then  the  optimum  weight  vectors  reduce  to 
some  quadratic  estimators  that  have  recently  been  derived  in  the  literature 
(ref. 6, 7).  In  these  cases  the  optimum  weighting  vectors  are  independent  of  the 
noise  crosspower  spectral  matrix  and  are  determined  by  the  array  geometry. 

This  work  is  part  of  a  continuing  RfiD  progranme  in  signal  processing  for  underwater 
acoustics  and  has  been  carried  out  under  tast  DST  79/069. 


2.  DERIVATION 


Let  X  =  {^jl 
be  w  =  {w^} 

is  defined*  by 


be  the  vector  of  outputs  of  the  K  receivers  at  some  frequency  f  and 
be  the  vector  of  receiver  weights.  The  covariance  matrix  R  (K  x  K) 


*  The  superscript 


denotes  the  Hermition  transpose  of  a  matrix  or  vector. 
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R  =  . 

In  the  absence  of  a  signal  the  receiver  noise  outputs  are  denoted  by  n  and  the 
noise  covariance  matrix  is  defined  by 

R  =  <  nn^  > 
n 


and 

R  =  Rc;  +  R  . 
S  n 


The  following  derivations  have  been  formulated  for  the  frequency  domain  but  can 
easily  be  extended  to  the  time  domain  by  a  suitable  redefinition  of  the  matrices 
and  vectors. 

Any  set  of  L  linear  constraints  on  the  weight  vector  w  can  be  expressed  as 

Mw  =  c  (1) 

where  M  is  the  L  x  K  constraint  matrix  and  c  is  the  L  x  1  vector  of  constraint 
values . 

2.1  Minimum  power 

This  criterion  chooses  the  weights,  w^,  such  that  the  total  output  power  ie 


w  Rw, 


is  minimised  subject  to  the  constraints  (1)  being  satisfied. 

Introducting  the  vector  X  of  undefined  Lagrangian  multipliers  the  cost 
function  H(w)  to  be  minimised  is  given  by 

HCw)  =  w*^  Rw  -  (Mw  -  c)^X  -  X^(Mw  -  c)  . 

PI 

Differentiating  the  above  equation  with  respect  to  w  and  equating  the  result 
to  zero*,  the  following  equation  is  obtained: 


Rw  -  M^X  =  0. 

(MPl 

Assuming  FT*  exists  then  the  optimum  weight  vector,  denoted  as  w^  ,  is 
given  by 


(2) 


*It  can  readily  be  proved  that  differentiating  with  respect  to  w  and  w  is 
equivalent  to  differentiating  with  respect  to  the  real  and  imaginery  components 
separately. 
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Furthermore  requiring  that  the  constraint  equation  (1)  is  satisfied  implies 
that 


(Mir'M^)A.  =  c 


(3) 


If  M  is  of  rank  L  (ie  of  full  rank)  then  it  follows  that 


X  =  (MJT  ‘  M^r  ‘  C 

and  hence 

(MP)  ^  IT  ‘  (MFT  '  M^)"  *  c 

which  has  been  derived  by  a  similar  argument  in  reference  2.  However,  more 
generally  it  follows  that  any  solution  to  equation  (3)  can  be  written  in  the 
form 


X  =  (Mir‘M^)'^c  +  (I  - 

where  z  is  an  arbitrary  vector,  and  denotes  the  Moore-Penrose  pseudoinverse 
of  A  (see  reference  8  or  Appendix  I  for  definition  of  A^) , 

In  this  case  equation  (2)  becomes 

^^CMP)  ^  ir*M^  (Mir‘M^)‘^c  +  IT'M^  (I  -  (MlT  '  M^) )  z 

which  reduces  to 


v^^(MP)  ^  R-iM^  (Mir*M^)‘"c  (4) 

since  (see  Appendix  I)  the  second  term  is  equal  to  zero.  The  output  power  is 
given  by 


(MP)* 


R  w 


(MP) 


(MR"  ‘  M^)  MR"  *  RR"  ‘  (MR"  ‘  M^)  ""c 
(MR"  ‘  M^ V  MR" '  M^‘  (MR-  ‘  M^)  *c 


=  c^(MR"'M^Vc 


by  virtue  of  equations  (I.l)  (ie  the  definition  of  the  Moore-Penrose 
pseudoinverse).  Note  that  both  the  expression  for  the  optimum  weighting 
vector  and  the  minimun  output  power  cai.  be  derived  directly  from  p.49  of 
reference  9. 
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2.2  Least  mean  square  derivation 

The  error  between  some  desired  output,  d,  and  the  weighted  receiver  outputs 
is  given  by 


6  =  d  -  W  X 


The  mean  square  value  of  this  error  can  be  shown  to  be  given  by 


H 


<  e^>  =  d“d  -  w”  <  d  X  >  -  <  d  >  w  +  w”  Rw 


H 


,  H  H  H  „ 

=  dd-w  v-v  w+w  Rw 


(5) 


where  the  covariance  between  d,  the  response  due  to  a  signal,  and  x,  the 
receiver  outputs  is  simply  the  steering  vector,  v,  if  the  signal  and  noise 
are  uncorrelated.  Thus  minimising  equation  (5)  subject  to  the  constraints 
implies  a  cost  function 


H(w)  =  dd  +  wRw-vw-wv  +  (c  -  Nrw)A+A  (c  -  Mw) 

Differentiating  with  respect  to  w^  and  equating  to  zero  implies  that  w 
is  given  by  ^ 

.  IT-  (lA  .  V)  ( 

Ensuring  that  satisfies  the  constraints  implies  that 


(MJT  ^ 


c  '  MlT  *  v 


and  hence  X  is  given  by 


X  =  (MIT^M^)'*'  (c  -  MirW)  +  (I  -  (Mir^M^))z 


where  z  is  arbitrary. 

Substituting  for  X  in  equation  (6)  it  can  be  shown  in  a  manner  similar  to 
that  used  in  Section  2.1  that 

^^(LMS)  ^  (Mir‘M^)''c  +  ir‘v  -  JT'm”  (MIT’M^)*  Mir‘v 

It  can  also  be  shown  that  the  output  power  is  given  by 

w  Rw  =  c”  (MFT'M^^J^c  +  ir‘v  -  IT'm'^  (MR" '  M^)  ""  MIT  > 

o  o 


v 
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2.3  Maximisation  of  beam  signal  to  noise  ratio 


If  the  signal  covariance  matrix  is  defined  by  R  then  the  beam  SNR  is  given 
by 


Sw 


H  D 

w  R  w 

s 

w  R  w 


Maximising  the  SNR  subject  to  constraint  (1)  implies  the  cost  function 


R  w  ,  H  H  . 

H(w)  » _ s  +  (c  -  w  M  jA  +  A  [c  -  Mwj 

H 

w  R  w 
n 

H 

Assuming  R^  to  be  non-singular,  differentiating  H(w)  with  respect  to  w  and 
equating  the  result  to  zero  implies  that 

P  w  .  -f(.  )  A 

o  o 


where 


P 


R  -  R  /g 
n  s  o 


and 


f(w) 


(w 


R  wP  /w^ 
n 


R  w 
s 


and 


„H(SNR)j^  ^  (SNR) 
o _ s  o _ 

^  „  (SNR) 

o  no 

By  an  analogous  derivation  to  the  previous  sections  it  follows  that 

v^CSNR)  ^  p- (Mr  *1^)" ‘c  (7) 

provided  both  P  and  Mr‘M^  are  non-singular. 
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Now  since  is  given  by  vv  then  it  follows  that 


vv 


P  =  R 


and  is  non-singular  unless  g^  -  v.  (From  p.29  of  reference  9 

I  P  I  =  I  R  I  (g  -  fT  *  v))  . 

n  ^o  n  ^ 

This  special  case  corresponds  to  the  single  constraint  of  unity  response  in 
the  look  direction  and  in  this  case  w  can  be  derived  by  the  method  in 

reference  1,  or  by  an  application  of  the  theorem  on  p.48  of  reference  9. 
Unfortunately  equation  (7)  expresses  w  as  a  function  of  g  which  in  turn 

( QMD  ^  ^  ^ 

is  defined  by  ,  In  Appendix  II  this  explicit  dependance  is  removed  and 

is  shown  to  be  given  by 


w 

o 


(SNR) 


where 


IT  '  (MIT  * 

n  ^  n 


c“(Mir^r>c 

-  _  ( IT  *  V  -  IT  *  NT  (MJT  ‘  NT)"  ' 

n  n  ^  n  ^ 

H 

c  y 


MR  v) 
n 


(8) 


y  =  (MJT  *  M^J  '  MR"  '  V 
^  ^  n  n 


As  is  also  shown  in  Appendix  II  che  optimum  gain  g^  is  given  by 


H  H 

H  c  y  y  c 

g  =  g  ~  y  MR*  '  V  +  „  „ 

c  (MR  M  )  c 


where  g 


is  defined  by 


fT  'v 


and,  as  discussed  earlier  (and  in 

optimum  processor  subject  only  to 
the  look  direction.  Furthermore 
one  of  the  L  multiple  constraints 


reference  1) ,  is  the  output  SNR  of  an 

the  single  constraint  of  unity  response  in 
if  unity  response  in  the  look  direction  is 
then  it  can  readily  be  shown  that 


\ 
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(MIT  *  M^*)"  ‘  c 

3.  RELATIONSHIPS  BETWEEN  THE  PROCESSORS 
All  optimum  weight  vectors  derived  are  of  the  form 

S'  *  M^(MS"  *  ‘  c  +  ocCS”  ‘  V  -  S"  '  M^(MS” ' m'”*)'  '  MS"  '  v) 

where 

(i)  S  =  R  and  oc  =  0  for  minimum  power  criterion, 

(ii)  S  =  R  and  a  =  1  for  least  mean  squares  criterion,  and 

c^(MR“^‘  M^)"‘c 

(iii)  S  =  R  and  a  =  for  maximum  signal  to  noise  ratio 

"  H 

c  y 

criterion. 

Now  if  it  can  be  shown  that 

-  5" 'm^3"  '  mS^'v  =  0  (9) 

then  it  follows  that  all  weight  vectors  are  of  the  form 


m“  (MS-‘ 


Expression  (9)  can  be  shown  to  be  equal  to  zero  provided  one  row  of  the 
constraint  matrix  M  is  equal  to  (or  a  multiple  of)  v^. 

This  can  be  simply  *  proved  by  noting  that  M  can  always  be  rearranged  such  that 
V  is  in  the  first  row.  Thus  by  taking  the  first  column  of  the  equation 


it  follows  that 


(MS”  ‘  *  (MS"  ‘  M^)  = 

(Ms”  ‘  '  (MS"  ‘  v) 


I. 


*The  author  is  indebted  to  Dr  A.K.  Steele  who  simplified  an  original  proof  of 
this . 
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Thus 


S'  '  (MS’  ’  M^f  ’  (M^"  *  v)  =  S'  ‘  /  1 

0 


=  S' '  V. 


Since  v  determines  the  look  direction  this  amounts  to  fixing  the  response  to 
known  plane  wave  signal. 

Thus  provided  the  same  constraint  matrix  M  is  used  to  estimate  w  and  w 

H  o  o 

and  one  row  of  this  matrix  is  given  by  v  then  it  follows  that 

w  (MP)  .  w 
o  o 

Furthermore  on  substituting 

R  =  R  +  vv^ 
n 

in  the  equation 

„  (MP)  =  JT  ’  (Mir  ‘  M^r  ‘  c 
o 


and  using  Woodbury’s 


Identity  for  (R^ 


H.-  1 
+  vv  ) 


it  follows  that 


(MP)  ^  IT  ‘  M^  (MIT*  M^f'c  +  a(ir‘  M^  (MR"’  M^f  '  MIT  ’  v  -  iT  *  v) 
o  n  n  n 


where  a  is  a  constant,  and  so 


(MP)  ^ 


n  '  n 


=  w 


(SNR)  . 


Thus  provided  the  same  constraints  are  imposed  and  a  fixed  look  direction 
constraint  is  imposed  it  follows  that  all  these  estimators  are  identical. 


a 

(IMS) 


4 .  EXAMPLES 


In  this  section  some  examples  of  the  processors  subject  to  some  obvious 
physical  constraints  are  given. 

4.1  Unity  (or  fixed)  response  in  the  look  direction 

It  can  readily  be  shown  that  for  a  processor  of  the  form  iT  * v  the  mean  beam 
power  output  decreases  as  the  signal  strength  increases .  This  can  be 
avoided  by  imposing  the  constraint  of  a  fixed  response  in  the  look  direction, 
ie 


H 

V  w 


1 


Thus 


M  = 


and 


c  =  1 


Then 


w  =  IT  ’  vCv^ir  'v)"  ’ 

o 

=  „  (LMS) 


since 


ir‘v  -  IT*  M^CMlT'M^r'Mir'v 


=  ir‘v  -  r“V 

=  0. 

Note  that  P  as  defined  in  Section  2.3  is  singular  and  hence 

w/snr)  .  r-ii'  A- ‘c 

is  ill -defined.  However  from  reference  1*  it  follows  that 


*Note  however  that  direct  substitution  in  equation  (8)  gives  the  correct 
solution  in  this  case. 


A 
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(SNR)  _  IT  'v 

(vV  ‘v)  . 

4.2  Fixed  response  in  the  look  direction  and  a  null  steered  in  a  specified 
direction 

A  particular  processor  that  has  received  some  attention(ref . 5)  maximises 
the  signal  to  noise  ratio  in  a  direction  v  while  steering  a  fixed  null  in 
another  direction  specified  by  u.  If,  as  in  the  previous  section  a  fixed 
(say  unity)  response  is  imposed  in  the  look  direction  then  the  constraint 
matrix  and  vector  are  given  by 


M  = 


and 


respectively. 

From  equation  (4)  it  follows  that 
W  =  R- >  (v,u) 


vU 


R"‘  (v,u)V‘  I  1 
0 


where 


=  (R-  *  v.R-  *  u)  /  SjjyN  '  '  /  1 


®vu 


g  =  V  IT  *  v,  g 

By  » 


R-'u. 


and 


Thus 


*  H  ^  1 

g  =  g  =  V  K  ‘  U. 
*UV  ^VU 


(MP)  ^  (frw,ir*u) 


^UV^VU^ 


g  -g  \  /I 
^U  ^uv 


^V 


L 


-  11  - 
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2  2  -  2  2 

®UV^VU 


Fr*v  -  g  ir'u 
*vu 


which  apart  from  the  scaling  factor  is  the  null  steering  weight  vector  derr 
by  Fenwick(ref . 5) . 


-Suv  \  /  1 


®u®v  ®uv®vu 


\  ®vu 


Thus  since 


it  follows  that 


^u^v  ^uv^vu 


^O 


2  2 
UV^VU 


which  again  agrees  with  the  expression  derived  in  reference  5. 
Also  it  can  readily  be  verified  that 


ir‘v  -  ir'M^  CMlT'M^r’  MIT* 


V  =  0 


thus  demonstrating  the  general  result  proven  in  Section  3  that 


^  (MP)  ^  ^  (LMS)  ^  ^  (SNR) 

o  o  o 


4.3  Fixed  response  in  look  direction  and  K-1  independent  nulls 

In  this  problem  all  the  degrees  of  freedom  are  incorporated  in  the  K  x  K 
non-singular  constraint  matrix,  M,  where  the  columns  of  M  are  the  phase  vcctv 
corresponding  to  the  K  directions. 
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It  then  follows  that 


^  (MR- > i/*)- ‘  c 

=  JT  ‘  *  RM"  ‘  c 

=  W'c  (10) 


Now  in  references  6  and  7  some  generalised  linear  processors  have  been  derived 
on  the  assumption  that  the  incident  distribution  consists  of  N  plane  waves. 
These  estimators  were  derived  on  the  basis  of  statistical  considerations  and 
deconvolution  of  the  array  response  respectively.  Defining  the  K  x  N  matrix 
whose  columns  are  the  phase  vectors  corresponding  to  the  N  assumed  directions 
as  V,  then  a  general  form(ref.6)  for  the  weighting  vector  in  the  i^^  direction 
is 


w . 
X 


=  V 


+H 


0 

1 

0 

6 


row 


/ 


(11) 


In  particular  when  N  =  K  and  the  directions  are  independent  then 


w. 

1 


-IH 


0 

i 

6 


.  th 

.^1  row 


which  as  can  readily  be  seen  is  a  special  case  of  equation  (10) .  The  fact  that 


H 

c 


(0, 


1, 


0) 


readily  demonstrates  the  observation  of  reference  6  that  for  a  given  steer 
direction  this  processor  steers  K-l  nulls  in  the  other  assumed  plane  wave 
arrival  directions . 

4.4  Overconstrained  case 

The  matrix  M  in  this  case  is  L  x  K  and  L>K.  To  obtain  consistent  solutions 
the  rank  of  M  must  be  less  than  or  equal  to  K.  If  M  is  chosen  to  be  of  full 
rank  then 


=  (M^)"  * 


and  hence 


13 
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m'^m  =  =  I 


The  L  X  L  matrix  MR“*m  is  singular  and  so  the  Moore-Penrose  pseudoinvcrsi' 
must  be  used.  Thus 


w 

o 


(MP) 


However  from  Appendix  I  it  follows  that 


(MIT  '  M^V  =  RM"" 

Thus 

»  .  r  ■  m"  m*"  rm*c 

o 

=  c 


=  (M^^r  '  Me 

As  in  the  previous  section  the  choice  of  M=V  and  c  =  (0,  1,  0) 

reduces  this  estimator  to  the  Least  Mean  Square  estimator  derived  in 
reference  6  and  to  a  special  case  of  the  estimator  derived  in  reference  7 
by  an  iterative  deconvolution  of  the  array  response. 


5 .  SUMMARY 

The  problem  of  constrained  optimisation  for  three  different  criteria  ie  Minimicv; 
Power,  Least  Mean  Squares  and  Maximum  Signal  to  Noise  ratio  has  been  solved  luui' 
some  very  general  conditions. 

In  all  cases  of  practical  interest  and  provided  identical  constraints  are 
imposed,  the  three  processors  have  been  shown  to  be  identical. 

A  number  of  examples  which  demonstrate  this  equivalence  have  also  been  given. 
Furthermore  these  examples  show  how  a  number  of  well  known  processors  are  all 
special  cases  of  the  constrained  optimum  processors.  In  particular  if  the 
weight  vectors  are  critically  or  over  constrained  then  weight  vectors  reduce  tn 
the  ideal  deconvolving  weights  which  have  been  derived  recently. 

This  paper,  in  common  with  previous  work  has  assumed  that  the  covariance  matrix 
is  non- singular.  This  is  not  always  necessary  and  further  work  is  proceeding 

to  derive  similar  expressions  which  do  not  require  non-singularity  of  the 
covariance  matrix. 
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APPHNDIX  I 

APPLICATION  OF  THE  MOORE-PENROSE  PSEUDOINVERSE 

1.1  Let  G  be  any  L  x  K  matrix  then  the  Moore-Penrose  pseudoinverse 
defined  by  the  following  equations: 

GG^G  =  G 

G^'GG'"  =  G* 

and  (GG"^)*^  =  GG* 

(G'^G)^  =  G*G 

Define  the  vector  u  by 

u  =  FT*  (1  -  (MIT'm^)'^  (Mir‘M^))z 

where  z  is  any  arbitrary  vector  and  R  is  a  K  x  K  positive  definite  mati  > 

and  M  is  any  L  x  K  matrix.  Denoting  (MIT'M^)  by  A  it  follows  that 

Ru  =  CI-A*A)^  MIT'riT’  m"  (I  -  A^\)  z 

=  (I-A'^A)^  ACI-A'^A)  z 

=  0  . 

Thus  since  R  is  positive  definite  it  follows  that 

u  =  0 

(Note  this  also  follows  directly  from  a  lemma  on  p.22  of  reference  10) . 

1.2  If  M  is  of  full  column  (ie  L  >  K)  rank  it  follows  that 

M'^M  =  =  I 

and  also  defining 

A  =  RM* 


WSRL-0178-TR 


16 


iv  follows  that 


M  =  i/*  = 


(I. 


and 


Thus 


and 


Also 


AA  =  MM"^ 


(M)“  =  M 


(AA)^  =  AA 


=  m'*'^  rm'^mm'^ 

=  A 


and  similarly 


AAA  =  A  , 


Hence  A  is  A  the  Moore  Penrose  pseudoinverse  of  A  where  A  =  (In 
this  case  it  can  readily  be  proved  by  direct  substitution  of  equation  (1.2) 
that 


(I  -  z  =  0 


=  0)  . 


and  hence  u 
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APPENDIX  II 

CONSTRAINED  MAXIMUM  SNR  PROCESSOR 


Since  P  is  defined  by 


where 


it  follows  that 


P  = 


R  -  R  /g 
n  s 


H  o 

w  R  w 
o  so 


w^  R  w 
o  no 


w^  P  w  =  0  , 

o  o 


Since  w^  is  given  by 


w  =  P"  *  1^  (MF  *  ‘  c 


this  then  implies  that 


(MF  ‘  M^r  '  c  =  0 . 


(II. 15 


If  R 


s 


is  equal  to  w 


H 


then  by  Woodbury’s  identity  it  follows  that 


where 


1 


F  *  +  F  ’  w“  F  ’ 
n  n  n 


8o  - 


(11.2) 


F  ‘  V . 


and 


WSRI.-0178-TR 


18 


Furthermore  it  holds  that 

(MF  ‘  *  =  (MIT^'  ‘  -  — 

P 

where 


y  =  (MIT  '  M^r  ‘  MR"  ‘  V 
^  ^  n  ■'  n 


and 


P  =  iT  ‘  (MlT  ’  '  MIT  *  V  +  a  . 

n  ^  n  ^  n 


Thus  equation  (II. 1)  becomes 


c“(Mir^*  m“)"  ‘  c 


H  H 
c  yy  c 


(v*^ir  ‘  (MIT  ‘  M^f  *  MIT  ‘  V  +  a) 
^  n  ^  n  ^  n  ^ 


Thus 


H  H 
c  yy  c 


a  = 


-  v^lT  ‘  (MR"  '  M^)"  ‘  MR"  ’  V . 
1  n  ^  n  n 


M*^)" 

Rearranging  equation  (II. 2)  and  substituting  the  above  expression  fora  gives 


g  =  g  -  v^R"  ‘  (MR"  ‘  M^r  ‘  MR"  *  V  + 


H  H 
c  yy  c 


^^o  '’v  n  '  n 


c^CMR" 


(II. 3) 


Hence 


w  =  r  *M^(Mr ‘M^r 'c 


r“  '  vv‘^r"  ‘ 

n  n 


yy 


|ir‘  +  ■ 

~1m”  j  (mr-'m”)-'  - - 1 

L" 

a 

J  L  "  P  \ 
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where  P  can  be  shown  to  be  given  by 


H  H 
c  yy  c 


c^CMlT^'  M“r‘c 


Thus 


w 

o 


IT  ‘  c 

n  ^  n  ^ 


fT^M^y 

n 


LS.  (y'^Mir^M^V  -  (?)ir’v 
pa  n  n 


=  IT ‘M^^CMIT  *1^)" 'c 
n  ^  n  ^ 

=  IT ‘m^CMIT ‘M^r ‘c 
n  n  ^ 


y"c  Cir‘v  -  ir‘M“  (MIT  ‘  M*^r  ‘  Mir  '  V  ) 

-  n  n  n  n 

c“(Mir  ‘  M^)"  ‘  c 

- - -  (IT  ‘  V  -  FT  *  Fr(Mir  *  Nr)“  ‘  Mir  ‘  v) . 

M  n  n^n-^n 
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